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Abstract. Wc analytically link three properties of nonlinear dynamical systems, 
namely sensitivity to initial conditions, entropy production, and escape rate, in z- 
logistic maps for both positive and zero Lyapunov exponents. We unify these relations 
at chaos, where the Lyapunov exponent is positive, and at its onset, where it vanishes. 
Our result unifies, in particular, two already known cases, namely (i) the standard 
entropy rate in the presence of escape, valid for exponential functionality rates with 
strong chaos, and (ii) the Pesin-like identity with no escape, valid for the power-law 
behavior present at points such as the Feigenbaum one. 
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1. Introduction 



At the onset of chaos, some traditional approaches to understand specific dynamical 
behavior, such as the entropy production, do not provide quantitatively nontrivial 
information on the state of the system. From the various physical phenomena which, 
at this limit, show a clear departure from the classical Boltzmann-Gibbs (BG) theory, 
the edge of chaos of one-dimensional maps has played the role of an archetypical system 
of study (see, for example, [U El |3] and references therein). In this manuscript we will 
focus on the relation between sensitivity to initial conditions, a possible escape within 
traps, and the loss of information in an one-dimensional map at the onset of chaos. 



1.1. The z-logistic map 

Let us focus on logistic-like maps x t+ i = f(x t ) (t = 0, 1, 2, ...) with 

f(x) = 1 - fx\x\ z (z > 1; < n < 2; - 1 < x < 1). (1) 

The expression of the sensitivity to initial conditions = limAx -s>o has been 
derived using exact analytic renormalization group theory (RGT), leading to 

6 = f = [1 + (1 - qsen)X qae j}^, (2) 

with q sen < 1, and e y = [1 + (1 — q)y] T ^ (e\ = e y ). When the Lyapunov exponent Ai is 
positive, we have q sen = 1. At the edge of chaos, from |4j, we have 

_ 1 (z- l)lna(z) 

q - 1- Qsen ~ In 2 ' [6) 
where a(z) is the ^-generalized Feigenbaum universal constant (see also [S]). In 
particular a(2) = 2.50290.... It follows that & = (1 + t) x <*<™ ~ t x ^ (t -»• oo). 
In this context, the nonadditive entropy [B] 

(Si = S B g = -^Pilnpi), 

i 

has been successfully used [7], where BG stands for Boltzmann-Gibbs, and ln q y = y 
(ln x y = lny). 



2. Pesin's theorem 



More recently, using RGT, it has been proved an extension of Pesin's theorem |Sj, valid 
for the entire family of these unimodal maps, making the important connection between 
the loss of information measure, via entropy production 
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and the sensitivity to initial conditions. This result can be expressed by the g-generalized 
Pesin-like identity 

with q ent = q sen . Whenever A x > we have 

Qent = Qsen = 1- (7) 

At the Feigenbaum point /i^ (where Ai = 0), Eq. ^ is satisfied. For z = 2 we have 
/ioo = 1.401155... (for z increasing from 1 to infinity, fioo{ z ) monotonically increases 
from 1 to 2). At this point, the pitchfork bifurcations accumulate, and the transition to 
chaos occurs. 

It is worth remembering that the Pesin equality relates the Lyapunov exponent 
Ai of nonlinear maps to the Kolmogorov-Sinai (KS) entropy /C. For relation (|6]), the 
probability distribution in Eq. Q is obtained as the frequency at which the i-th cell in 
the state space, i = 1,2, is occupied at a given time. The main difference that one 
would expect between a possible g-generalization of the Kolmogorov-Sinai rate, noted 
fC q , and the entropy production rate K q , in Eq. (|6]), is that in the latter the initial 
conditions for the trajectories have a specific distribution Pi(0) while the g-KS-entropy 
case, fC q , considers all the trajectories from their initial positions to the time t — > oo. 
The connection between K\ and /C has been studied for several chaotic maps [9] and it 
has been determined that typically the asymptotic equality Ki = K, holds. In numerical 
calculations, this corresponds to an intermediate stage of evolution, after a transient and 
before reaching the asymptotic (saturated) state. In the same manner, we are interested 
here in the temporal region where K q = K q , and, as shown in [8], it will be necessary 
to study only uniform initial distribution. 

We introduce the definition of the sensitivity factor £ t for each time subsequence 
k, k = 1,2,3..., and for two trajectories x tk and y tk with initial condition Xi n , yi n 

it = lim \ytkiVin) — x t k {Xin)\ 

lifc„-xi„|->o \y tk =o{yin) - x tk =o{x in )\' 
We will consider in this map a partition of W equal-length boxes of size A, and an 
ensemble of iV trajectories with initial condition near Xi n = 1 (in what follows, we will 
omit the k- index, that, as we mentioned, describes a possible time subsequence). Then, 
the entropy production can be calculated considering the occupancy frequency of each 
box. The total number of cells that the ensembles of trajectories occupy at time t will be 
W t = A f /A, where A t is the total size occupied by the ensemble. As pointed out in [8], 
in the limit of A — > 0, the total number of occupied cells will be uniformly distributed 
as Wt = C,t- This result implies that the total entropy at time t will be simply given by 
S q = \ia q W t = X q t, leading to Eq. (6). In our analysis we will consider the same type 
of expansion for a generic dynamic map, that is 

W t = ^ (8) 

This relation is based on maps that have been studied in the literature, as discussed 
above, some of them homomorphism of important continuous chaotic dynamical systems 
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[TO] . The key point in our discussion will be the relation between this expansion and 
the functionality of the escape rate. 

3. Open systems 

In the presence of escape gates in the state space, the exponential decay behavior 
(corresponding to Ai > 0) for a given number of initial trajectories in the phase space has 
been extensively reported [JJ. This means that for a given number of initial trajectories, 
the number of them that remain in the system at time t, N t , will be given by 

N t = N e^ *. (9) 

For a dynamical system that has an exponential uniform expansion, and under the 
presence of an escape rate that behaves as the previously mentioned one, the number 
of occupied cells at time t can be expressed as 

W t = e" 71 *e Al *. (10) 
And using Eq. Q, with q = 1, we can obtain the well known relation 

S 1 = hy[W t ] = (\ 1 -j 1 )t. (11) 
Consequently, given that K 1 = lim^oo sx, we have 

#i = Ai-7i, (12) 

which connects, for such systems, the entropy production K x , the Lyapunov exponent 
Ai, and the escape rate 71. 

Notice that in principle the functionality of both the divergence of initial conditions 
and the escape should be the same in order for this relation to exist (exponential in the 
strong chaos, power-law at its onset). 

Consider now an uniform expansion that behaves as a power law (instead of as an 
exponential), for example, the logistic map at ix^: what will the entropy production be 
under a power law-type of escape? Note that, if the decreasing number of trajectories 
in the presence of escape is exponential, the behavior of the occupied bins after some 
time will be far from power law. 

Assuming that our system is such that the escape asymptotically follows 

N t oc £ _7W , (13) 
and taking into account that 

N t = N e~2 q c eact oc £~V(9«>c-i) (t->oo), (14) 
we can identify 

1 . . 

Iiesc = ~ — T7 • (15) 

Hesc J- 

The production of occupied bins, W t , and sensibility £ t (see Eq. (J8l)) will behave as 
W t = & oct^t -7 *". (16) 
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Then, finally, the entropy production (from Eq. (15])) is given by 
with 

K. 



Qent -i 

J- ~~ Qent 



Eq. (17) generalizes Eq. (12), the well known relation valid for chaotic systems (i.e., 



Ai > 0) in the presence of escape (see for instance pQ). Notice that the q indices present 



in Eq. (17) do not need to be equal among them. Indeed the exponents in the power 



laws shown in Eqs. ^ and (13) do not necessarily coincide. This implies, through Eq. 



(18), that also q en t generically differs from them. 



3.1. Example 

As a concrete example of a system that has the properties discussed above, we can 
mention the following non-autonomous map 

x t +i = [1 - ^(x 2 t + y 2 )] cos e t 

y t +i = [1 - fi(x 2 t +y?)}sme t , (19) 

with 

In (t + 1) 

which is effectively a rotating logistic map. The values of cosBj (sinBj) will be ±1 
(0) exactly when the map expands in the x axis (with y = 0). This system can 
simultaneously show a critical behavior such as that at the onset of chaos (i.e., at 
/Zoo), and power-law escaping due to holes that are uniformly distributed on the line 
y = 0. In such situation we can achieve an exact solution for the entropy production, 
using the usual ensemble of initial trajectories at y — 0, near x = 1 |8], with t = 2 l — 1, 
i = 1,2,3.... 

Then, if we introduce a fraction of (1 — 5) traps (or holes) in the domain 
— 1 < x < +1 , equidistantly distributed in the line y — 0, the system will suffer a 
decrease in the number of trajectories in a proportion given by 5 l . The number of 
occupied bins, at time t = 2 l — 1, will be given by 

W(t) = W(2 i - 1) = (1 + t)Ven +£-§. (21) 

Finally, after a straightforward calculation of the entropy production with 

1 

qent — 1 - ~ ln(5) > ( 22 ) 

A qsen + ln(2) 

we obtain a generalization of the Pesin relation of rate parameters, including power-law 



escaping Eq. (17) where we have used Eq. (18), and 



In 5 

7 9esc = -^. (23) 



In Fig. 1 we illustrate Eq. (13), and in Fig. 2 we show a corroboration of equality (17). 



Sensitivity, entropy, and escape rates at the onset of chaos 



6 




Figure 1. Fraction of points, n(t) = Ng/Nt, remaining in the system versus time 
using 5 — 6/7 ~ 0.86 and z = 2, in log-log scale. Nq = 10 6 uniformly taken within the 
interval [1 — 10~ 10 , 1]. The fit, dashed line, shows a escape parameter 7 9esc = 0.216... 
while the theoretical one, calculated from Eq. (23), is ~f qesc = 0.2223.... 




70 S q (t) 140 



Figure 2. Sensitivity to initial condition versus entropy production, see Eqs. (16) 



and (17), for different values of z. For z = 2 and 5 = 0: 
lent = Qsen = 0.244...; while for z = 2 and 5 = 6/7: 7 g 
= 1.1012 



K„ 



ni - \ 

0.222.. 



n = !- 32 - 

from Eq. 



and 



K, 



and q e 



0.0919..., from Eqs. (17) and (18). Similar results can 



(23) 



be obtained for the other values of z. The holes are uniformly distributed in the line 
y = 0. The continuous line correspond to a fit with a slope 1.004..., numerically very 
close to unity, as expected. These examples neatly illustrate the validity of Eq. (17): 
the ordinate corresponds to (A gsE?l — 7 9eac ) t, and the abscissa corresponds to K qent t. 



4. Final remarks 



To conclude, let us mention that one of the important aspects of dynamical systems 
with escape is their connection with physical scenarios that present boundaries in its 
phase space [11]. Our results show that, at the edge of chaos, a thermodynamical- 
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like quantity (entropy production) and a dynamical quantity (sensitivity to the initial 
conditions) are connected even in the presence of escape. This type of connection, well 
known for chaotic dynamical systems, is then valid, and natural, when considering a 
suitable entropy that takes into account the particular regime (multifractal structure in 
phase space, power-law escaping, etc) under which the system evolves. 

The logistic map dynamics at the Feigenbaum point has been recently studied in 
much detail [8J. Its structured behavior makes possible an analytic treatment. The 
results presented here used this amenable characteristic of this particular map, but our 
main point goes beyond it. Our results would be applicable in principle to any system 
where both trajectory divergence within the phase space, and escaping to outside of the 
system, follow power laws. Typical examples would be deterministic anomalous diffusion 
in Lorenz gas, in Hamiltonian systems (e.g., those numerically explored in [12]), and 
spatio-temporal intermittency in dissipative systems. As pointed out, our example, 
at the onset of chaos, can be an archetypical illustration where the phenomenology 
discussed here occurs, but also more complicated systems where attractors (see also [13] ) 
with low-dimensional behavior and its transients, that connect them, are embedded in 
a larger system. 
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